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Abstract 

\^ ' We present results from a study of subtractive renormalization of the 

00 ■ fermion propagator Dyson-Schwinger equation (DSE) in massive strong- 

(N 
O 



X 



couphng quenched QED4. Results are compared for three different fermion- 
photon proper vertex Ansdtze: bare j^, minimal Ball-Chiu, and Curtis- 
Pennington. The procedure is straightforward to implement and numerically 
stable. This is the first study in which this technique is used and it should 



Q_i' prove useful in future DSE studies, whenever renormalization is required in 

^ • numerical work. 



I. INTRODUCTION 



Strong coupling QED in three space and one time dimension has been studied within the 
5^ ! Dyson-Schwinger Equation (DSE) formalism for some time 0-0, in order to see whether 
there may be a phase transition to a nontrivial "local" theory at high momenta as a 

model for dynamical chiral symmetry breaking (DCSB) in walking technicolor theories 
and also as an abelianized model for nonperturbative phenomena in QCD [§,0. For a recent 
review of Dyson-Schwinger equations and their application see Ref. [|ll|]. The usual approach 
is to write the DSE for the fermion propagator or self-energy, possibly including equations 
for the photon vacuum polarization |12,13] or the fermion-photon proper vertex [Q]. An 



Ansatz is made for the undefined Green's functions that contain the infinite continuation 
of the tower of DSE's. The resulting nonlinear integral equations are solved numerically in 
Euclidean space |T^ by iteration. DCSB occurs when the fermion propagator develops a 
nonzero scalar self-energy in the absence of an explicit chiral symmetry breaking (ECSB) 
fermion mass. We refer to coupling constants strong enough to induce DCSB as supercritical 
and those weaker are called subcritical. We write the fermion propagator as 

S(v) = ^(^'^ = 1 fl) 
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with Z{p^) the finite momentum-dependent fermion renormahzation, and B{p^) the scalar 
self-energy. In the massless theory (i.e., in the absence of an ECSB bare electron mass) 
by definition DCSB occurs when B^p^) ^ 0. Note that A^p^) = l/Z{p^) and M^p^) = 

Many studies, even until quite recently, have used the bare vertex as an Ansatz for the 
one-particle irreducible (1-PI) vertex r^(fc,p) ||,|,|gJ|JT|JT|JTl , despite the fact that this 
violates the Ward-Takahashi Identity (WTI) []TB|. It is also common, especially in studies 
motivated by walking technicolor theories 0, to find vertex Ansdtze which claim to solve the 
WTI, but which still possess kinematic singularities in the limit of zero photon momentum 

= {k — pY —>■ [|1^,|T^. With any of these Ansdtze the resulting fermion propagator 
is not gauge-covariant, i.e., physical quantities such as the critical coupling for dynamical 
symmetry breaking, or the mass itself, are gauge- dependent |T^. A general form for r'^(A;,p) 
which does satisfy the Ward Identity was given by Ball and Chiu in 1980 it consists 



of a minimal longitudinally constrained term which satisfies the WTI, and a set of tensors 
spanning the subspace transverse to the photon momentum q. 

Although the WTI is necessary for gauge-invariance, it is not a sufficient condition; fur- 
ther, with many of these vertex Ansdtze the fermion propagator DSE is not multiplicatively 
renormalizable, which is equivalent to saying that overlapping logarithms are present. There 
has been much recent research on the use of the transverse parts of the vertex to ensure 
both gauge-covariant and multiplicatively renormalizable solutions []20|-p^, some of which 
will be discussed below. 

What is common to essentially all of the studies so far is that the fermion propagator is 
not in practice subtractively renormalized. Most of these studies have assumed an initially 
massless theory and have renormalized at the ultraviolet cutoff of the integrations, taking 
Zi = Z2 = 1. Where a nonzero bare mass has been used PJT^, it has simply been added 
to the scalar term in the propagator. Although there have been formal discussions of the 
renormahzation P,|5| J23|J2B|] , the important step of subtractive renormahzation has not been 
performed. 

We describe here the results of a study of subtractive renormahzation in the fermion 
DSE, in quenched strong-coupling QED4. (In the context of this study of QED, the term 
"quenched" means that the bare photon propagator is used in the fermion self-energy DSE, 
so that Z3 = 1. Virtual fermion loops may still be present, however, within the vertex 
corrections.) We believe this is the first calculation in which the subtraction is performed 
numerically as the integral equations are being iterated. This new technique should prove 
useful in other DSE studies. Results are obtained for DSE with three different vertices: the 
bare 7^^, the minimal Ball-Chiu vertex form |T^], and the Curtis- Pennington vertex 

Organization of the paper is as follows. The DSE for the renormalized fermion prop- 
agator, and Ansdtze for the proper vertex, are discussed in section |T|. The subtractive 
renormahzation is described in section |ITT[ Section [IV A| presents results at subcritical cou- 
plings, while Section [IV B| discusses the results for couplings above the DCSB threshold. 
Conclusions and discussion are given in section V. 



II. DSE AND VERTEX ANSATZE 



The DSE for the renormalized fermion propagator, in a general covariant gauge, is 
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S-\p') = Z,{f,,A)[^-mo{A)]-tZ,{f^,A)e'f-0^^rS{k)T''{k,p)^^^^ (2) 

here q = k —p is the photon momentum, n is the renormahzation point, and A is a regular- 
izing parameter (taken here to be an ultraviolet momentum cutoff). We write mo (A) for the 
regularization-parameter dependent bare mass. In the massless theory (i.e., in the absence 
of an ECSB) the bare mass is zero, mo (A) = 0. The physical charge is e (as opposed to the 
bare charge eo), and the general form for the photon propagator is 

with ^ the covariant gauge parameter. Since we will work in the quenched approximation 
and the Landau gauge we have = Cq = 47rQ;o and 

for the photon propagator. 



A. Vertex Ansatz 

The requirement of gauge invariance in QED leads to the Ward-Takahashi Identities 
(WTI); the WTI for the fermion-photon vertex is 

q,T^{k,p) = S-\k)-S-\p) , (5) 



where q = k —p p9| , |3Q| . This is a generalization of the original differential Ward identity, 
which expresses the effect of inserting a zero-momentum photon vertex into the fermion 
propagator, 

dS-\p) 

^^ = r(p,p). (6) 

In particular, it guarantees the equality of the propagator and vertex renormahzation con- 
stants, Z2 = Zi. The Ward-Takahashi Identity is easily shown to be satisfied order- by-order 
in perturbation theory and can also be derived nonperturbatively. 

As discussed in |Tl| , |31| , this can be thought of as just one of a set of six general require- 
ments on the vertex: (i) the vertex must satisfy the WTI; (ii) it should contain no kinematic 
singularities; (iii) it should transform under charge conjugation (C), parity inversion (P), 
and time reversal (T) in the same way as the bare vertex, e.g., 

C-'T,{k,p)C = -Tl{-p,-k) (7) 

(where the superscript T indicates the transpose); (iv) it should reduce to the bare vertex in 
the weak-coupling limit; (v) it should ensure multiplicative renormalizability of the DSE in 
Eq. (^; (vi) the transverse part of the vertex should be specified to ensure gauge-covariance 
of the DSE. 
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Ball and Chiu [T^ have given a description of the most general fermion-photon vertex 
that satisfies the WTI; it consists of a longitudinally-constrained (i.e., "Ball-Chiu") part 
Tf^Q, which is a minimal solution of the WTI, and a basis set of eight transverse vectors 
T/'(/c,p), which span the hyperplane specified by g^Tf (/c,p) = 0, q = k — p. The minimal 
longitudinally constrained part of the vertex is given by 



jk + pY 

^2 _ p2 



+ 



- [Bik') - B{p')] 



The transverse tensors are given by 



Trik,p 
mk,p 

T,'{k,p 
T^{k,p 

T^{k,p 
Ts'iKp 



P^{k ■ q) - k^{p ■ q) , 

q'Y-q'^, 
p'k^a^xT,^ , 

l'{k'-p')-{k + pr ^, 
{k^-p^^ 



b^'i^ +^)-{k + pY] + {k + pfp^k^a^x 



--f^'p^k'^axu -p^ ^-k^" 
A general vertex is then written as 



(8) 



(9) 
(10) 

(11) 
(12) 
(13) 
(14) 

(15) 
(16) 



r^(A;, p) = r^^( A;, p) + J] r, (P, p^ q')Ttik p) 



i=l 



(17) 



where the Tj are functions which must be chosen to give the correct C, P, and T invariance 
properties. 



The work of Curtis and Pennington p3|-p6[] was mentioned above in connection with 
the specification of the transverse vertex terms in order to produce gauge-invariant and 
multiplicatively renormalizable solutions to the DSE. In the framework of massless QED4, 
they eliminate four of the transverse vectors since they are Dirac-even and must generate 
a scalar term. By requiring that the vertex T^{k^p) reduce to the leading log result for 
k ^ p they are led to eliminate all the transverse basis vectors except T^, with a dynamic 
coefficient chosen to make the DSE multiplicatively renormalizable. This coefficient has the 
form 



Teie,p^q') = ^[Aie) - Aip')]/dik,p) 



(18) 



where d{k,p) is a symmetric, singularity-free function of k and p, with the limiting behavior 
\imk2 d{k,p) = k"^. [Here, A{p^) = l/ZijP') is their l/jF(p^).] For purely massless QED, 
they find a suitable form, d{k,p) = [k"^ —p^Y I ^ This is generalized to the case with 

a dynamical mass 

M(p2), to give 
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d{k,p) 



^2 _j_ p2 



(19) 



They establish that multiphcative renormahzabihty is retained up to next-to-leading-log or- 
der in the DCSB case. Subsequent papers establish the form of the solutions for the renor 

751 and studied the gauge-dependence of the solutions [56 



Dong, 



malization and the mass 
Munczek and Roberts subsequently showed that the lack of exact gauge- covariance of 
the solutions is in part due to the use of a hard cutoff in the integral equations. The fact 
that there is still some residual gauge-dependence in the physical observables such as the 
chiral critical point shows that the C-P vertex Ansatz is not yet the ideal choice. Dong, 



Munczek and Roberts [271 derived an Ansatz for the transverse vertex terms which satisfies 



the WTI and makes the fermion propagator gauge- covariant. 

Bashir and Pennington [^] have recently described a vertex Ansatz which makes the 
fermion self-energy exactly gauge- covariant, in the sense that the critical point for the chiral 
phase transition is independent of gauge. Their work is a continuation of that of Dong, 
Munczek and Roberts, and indeed their vertex Ansatz corresponds to the general form sug- 
gested in [^ . However, the kinematic factors T2^3,6,8 are extremely complicated in form and 



depend upon a pair of as yet undetermined functions Wi^2{k'^ , P^) which must be chosen to 
guarantee that the weak-coupling limit of P'^ matches the perturbative result. Renormal- 
ization studies of the DSE using this vertex Ansatz should be interesting and represent a 
direction for further research. 

In any case, for the solutions to the fermion DSE using the C-P vertex, the critical point 
for the chiral phase transition has been shown to have a much weaker gauge- dependence 
than that for the DSE with the bare or minimal Ball-Chiu vertices [|3^. In this work we will 



primarily use the Curtis-Pennington Ansatz but for the sake of comparison will also present 
results for the bare vertex. Some solutions are also obtained with the minimal Ball-Chiu 
vertex which, like the Curtis-Pennington vertex, satisfies the WTI, but which does not lead 
to approximately gauge-invariant solutions nor multiplicative renormahzabihty. 

For each vertex Ansatz, the equations are separated into a Dirac-odd part describing the 
finite propagator renormalization A{p'^), and a Dirac-even part for the scalar self-energy, by 
taking |Tr of the DSE multiplied by and 1, respectively. The equations are solved 

in Euclidean space and so the volume integrals / d'^k can be separated into angle integrals 
and an integral / dk"^; the angle integrals are easy to perform analytically, yielding the two 
equations which which will be solved numerically. 



III. THE SUBTRACTIVE RENORMALIZATION 

The subtractive renormalization of the fermion propagator DSE proceeds similarly to 
the one-loop renormalization of the propagator in QED. (This is discussed in [|ll| and in 



33| , p. 425ff.) One first determines a finite, regularized self-energy, which depends on both 



a regularization parameter and the renormalization point; then one performs a subtraction 
at the renormalization point, in order to define the renormalization parameters Zi, Z2, 
which give the full (renormalized) theory in terms of the regularized calculation. 

A review of the literature of DSEs in QED shows, however, that this step is never 
actually performed. Curtis and Pennington [EB] for example, define their renormalization 
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point at the UV cutoff. Miransky ||^ gives a formal discussion of the variation of the mass 
renormahzation ^^(/i, A), but does not implement it numerically. 

Many studies take Zi = Z2 = 1 [ pO| , p^ -p6|; this is a reasonable approximation in cases 
where the coupling is small (i.e., ~ 1-15), but if Oq is chosen large enough, the value of 
the dynamical mass at the renormahzation point may be significantly large compared with 
its maximum in the infrared. For instance, in [^], figures for the fermion mass are given 
with ao = 0.97, 1.00, 1.15 and 2.00 in various gauges. For ao = 2.00, the fermion mass at 
the cutoff is down by only an order of magnitude from its limiting value in the infrared. 

Repeating the arguments from |]TT|, one defines a regularized self-energy A; p), lead- 



ing to the DSE for the renormalized fermion propagator. 



S-\p) = Z,ifi,A)y-mo{A)] - A;p) 

= ^-m(/i)-S(/i;p), (20) 



where the (regularized) self-energy is 

S'(/x, A; p) = zZi(/i, A)e' f -^J-^Sifi; fc)r (/i; k, p)D,,{fi; {p - k)) . (21) 

[To avoid confusion we will follow Ref. |TT| and in this section only we will denote regularized 
quantities with a prime and renormalized ones with a tilde, e.g. S'(yU, A; p) is the regularized 
self-energy depending on both the renormahzation point fi and regularization parameter 
A and Ti{fi]p) is the renormalized self-energy.] As suggested by the notation (i.e., the 
omission of the A-dependence) renormalized quantities must become independent of the 
regularization-parameter as the regularization is removed (i.e., as A ^ 00). The self-energies 
are decomposed into Dirac and scalar parts, 

S'(/x, A; p) = S:,(/x, A; p^) ^ + S;(/i, A; p') 

(and similarly for the renormalized quantity, S(/i,p)). By imposing the renormahzation 
boundary condition. 



one gets the relations 



for the self-energy. 



S''iP)^,_.,=I^-m{f^), (22) 



S,,,(/i;p^) = S:,,,(/x, A;p^) - S:,,,(/i, A; /x^) (23) 



Z2(/i,A) = l + S:,(/i,A;/i2) (24) 

for the renormahzation, and 

mo(A) = [m(/i) - S;(/i, A; /i^)] /Z^ifi, A) (25) 

for the bare mass. In order to reproduce the case with no ECSB mass, for a given cut- 
off A, one chooses m(/i) = S^(yU, A;/i^) so that the bare mass mo(A) is zero. The mass 
renormahzation constant is given by 
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Z„(/i,A) = mo(A)/m(/i), (26) 

i.e., as the ratio of the bare to renormahzed mass. 

The vertex renormahzation, Zi{fi, A) is identical to Z2{fi, A) as long as the vertex Ansatz 
satisfies the Ward Identity; this is how it is recovered for multiplication into S'(;U, A;p) in 
Eq. (pri). It will be noticed that this is inappropriate for the bare- vertex Ansatz since it fails 
to satisfy the WTI; nonetheless, since for the bare vertex case there is no way to determine 
Zi(/i, A) independently we will use Zi = Z^ for the sake of comparison. In the Landau gauge 
for the bare vertex these will then both be 1, since in this case A;p^) = for all as 



is well known [11 



We carry out this step of the renormahzation, using the bare vertex 7^^, the minimal 
Ball-Chiu vertex of Eq. (P), and that with the Curtis- Pennington transverse term added [see 
Eqs. (|TB|,|IUp]. The equations are again first written in Minkowski space, and then rotated 
to Euclidean space. These are then solved by iteration from an initial guess, with a wide 
range of cutoffs, and with various renormahzation points /i^ and renormahzed masses m(/i). 



IV. RESULTS 

Solutions were obtained for the DSE with the Curtis-Pennington and bare vertices, for 
couplings «o from 0.1 to 1.75; solutions were also obtained for the minimal Ball-Chiu vertex, 
with couplings from 0.1 to 0.6 (for larger couplings the DSE with this vertex was suscepti- 
ble to numerical noise). In Landau gauge, the critical coupling for the DSE with bare vertex 
is a^^''*' = 7r/3; the critical coupling for the Curtis-Pennington vertex is = 0.933667 



3^ , and that for the Ball-Chiu vertex is expected to be close to these two values. Since 
there are some qualitative differences between the solutions with subcritical and supercritical 
couplings, we shall describe these separately. 



A. Subcritical Couplings 

Solutions for the Curtis-Pennington and bare vertex Ansdtze were obtained for values of 
the coupling ao from 0.1 to 0.9. The program using the Ball-Chiu vertex was also run, with 
couplings ao from 0.1 to 0.6. For the solutions in the subcritical range, the renormahzation 
point fi'^ = 100, and renormahzed masses were either m{fi) = 10 or 30. Ultraviolet cutoffs 
were 1.0 x 10^^ ^nd 1.0 x 10^^ 

The family of solutions for the Curtis-Pennington vertex with m(/i) = 10 is shown in 
Fig. |l]. Fig. H a) shows the finite renormahzation A(p^); note that for the solutions with 
bare 7^^, we would have A(p^) = 1 for all values of the coupling. Fig. |l]b) shows the masses 
M(p2). 

Since the equations have no inherent mass-scale, the cutoff A, renormahzation point n, 
m{n), and units of M{p'^) or B{p'^) all scale multiplicatively, and the units are arbitrary. 
For instance, the solutions for any given coupling in Fig. |I] could represent a particle with 
mass lOMeV at the renormahzation point jj? = — lOOMeV^, with the units of in MeV^ 
and those of B in MeV; or it could represent one with units in GeV, or in electron masses. 
Since in four dimensions the coupling has no mass dimension, it would remain unchanged 
for all such choices of units. This is true for all the graphs that follow. 
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A comparison of the solutions for the three different vertices, for the couphng ao = 0.5, 
IS given m Fig. |. As can be seen in Fig. ^ a), the finite renormahzation A{p'^) deviates 
less from unity in the Curtis- Pennington case. The scalar self-energy curves B{p^), shown 
in Fig. H b), look qualitatively similar. However, since M{p'^) = B{p'^)/A{p'^), there are 
differences both in the infrared value of mass (for which we have Mcp < Mbc < ^bare for 
the same ao), and in the asymptotic scaling. 

The slight difference in scaling of the mass functions for different vertex Ansdtze is shown 
in Fig. 1^, for = 0.5. Miransky has shown ^ that in Landau gauge with the rainbow 
approximation, the fermion mass scales asymptotically as M(p'^) ~ {p'^)~^^~^'^ with the 
anomalous dimension 7'(ao) = — ao/i^Tr/S). For ao = 0.5, the anomalous dimension is 
7' = .3614329, giving the power | — 7' = 0.138567. When all the solutions are multiplied 
by (p2)5-7'(o.5)^ difference in anomalous dimensions shows up as a difference in slope 
of the mass curves on a log- log plot. The anomalous dimensions for the mass functions 
with Ball-Chiu and Curtis-Pennington vertices are found graphically, to be approximately 
7bc = 0.35836 giving the power | — 7cp = 0.14164, and 7cp = 0.35843 giving the power 
I — 7cp = 0.14157, respectively. The small dips apparent at the end of two of these curves 
are due to having a hard cut-off A^ = 10^^. As A is increased these move to higher momenta 
also. 

For the equations with the Curtis-Pennington vertex, the renormahzation constants Z2 
and Zm, and the bare mass mo(A) were evaluated for = 0.1 through 0.9, using cutoffs 
1 X 10^^ and 1 x 10^^ and renormahzation conditions /x^ = 100, m{n) = 10. These are 
presented in Table |. Specifically for ao = 0.5, the renormahzation constants Z2(/i, A) and 
Zm{fJ',A.) were calculated for several values of /i, using cutoffs A^ = 1 x 10^^ and 1 x 10^*. 
These are given in Table 0, and the mass renormalizations Zm are plotted as a function of 
/i^/A^ in Fig. |^. All these solutions used bare masses chosen so as to duplicate the mass 
solution for /i^ = 100, m(/i) = 10. It can be seen that as the value of /i increases, Z2 
approaches 1, and also that Z2 is almost independent of the cutoff A, at least for A ^ /i. 
The mass renormahzation also scales asymptotically (i.e., for sufficiently large A and fi) as 
ZmifJ', A) oc (yU^/A^)2~'''cp. The slight deviation from this scaling behaviour evident in Fig. ^ 
occurs in the region where /i is not sufficiently large. 

The solutions are extremely stable as the cutoff is varied. The ultraviolet momentum cut- 
off is varied by six orders of magnitude with the mass unchanged within numerical accuracy 
for sufficiently large A. 

B. Results above DCSB Threshold 

The program using the Curtis-Pennington vertex was run with several values of the 
coupling from 0.97 to 1.75, and with several different cutoffs. In all cases the mass and 
finite renormahzation were stable with respect to very large variations in cutoff. Typical 
results for the self-energy B{p'^) are given in Fig. ^. These were generated using = 1.00, 
/i^ = 100, and m(/i) = 10 (in arbitrary units). In both graphs, the cutoff is stepped by 
factors of 100 over eight orders of magnitude, and again the variation in the solutions is of 
the order of the tolerance of the calculation. 

The bare- vertex DSE was run with couplings ao = 1.15, 1.5 and 1.75 for comparison 
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to the Curtis-Pennington solutions. A comparison of typical mass curves is given in Fig. |^. 
Here the mass functions are compared directly since for the bare vertex the scalar self-energy 
is exactly the mass. 

An unexpected feature of the equations (for supercritical couplings only) is that for any 
set {ao, /i^, m(/i)}, as the cutoff is increased there is a small region where the dynamical mass 
becomes negative. This behavior was verified to be insensitive to number of grid points, and 
was stable as the cutoff was increased above the region of the negative self-mass, so it seems 
to be a real effect. The significance of this for QED is not completely understood. One 
possibility is that it may signal the the failure of multiplicative renormalizability for the 
model DSE, in which case further refinements of the vertex Ansatz may be called for. It 
should be emphasized that this negative dip in the scalar self-energy is very small unless the 
coupling ao approaches 2. For instance, for ao = 1-00, the negative peak was ~ 2.4 x 10~^ 
the size of the renormalized mass. Specifically, we found that for /i^ = 100, m(/i) = 10, 
the values are S(0) = 10.745, 5(1.634 x 10^^) ^ -.00024; for the case with = 1-75, 
/i2 = 1 X 10^ m(/i) = 500, the values are 5(0) = 531.8, 5(3.338 x 10^) = -28.3. An 
example is shown for ao = 1.75, /i^ = 1 x 10^, m(/i) = 500 in Fig. 0. 

There is no apparent scaling behavior for the mass renormalization Zm in the super- 
critical case. The renormalization constant Z2{fi,A) remains finite and well-behaved with 
increasing A as expected in the Landau gauge. Typical values, using ao = 1.15 and the 
Curtis-Pennington vertex, are given in Table |T|. From the results in this table and from 
Eqs. ( p5| , p6D we see that Sg(/i, A; /i) first increases beyond m(/i) as A increases so that mo(A) 
goes negative. As A is further increased the last line of this table indicates that S'^ then 
begins to decrease. 



V. SUMMARY AND CONCLUSIONS 

We have described preliminary results in a study of four- dimensional quenched QED, 
with subtractive renormalization performed numerically, during the calculation. We believe 
that this is the first calculation of its kind, and the technique described here will be appli- 
cable elsewhere (e.g., also in QCD), whenever numerical renormalization is required. The 
importance of this approach is that it allows studies using different regularization procedures 
(e.g., various soft versus hard momentum cut-offs) and/or different cut-offs to be meaning- 
fully compared at the same renormalization point. This includes also comparisons with 
results from lattice studies of QED, which should prove useful in providing further guid- 
ance in the choice of reasonable Ansdtze for the vertex and photon propagator. Without 
renormalization only the unrenormalized, regulated quantities would be be obtained and any 
such comparisons would be meaningless. In addition, in order to study the nonperturbative 
behaviour of renormalization constants such as Zi(/i, A), Z2(/i, A), and Zm(/U, A) they must 
be numerically extracted and so a method such at that described here would be essential. 

The Curtis-Pennington vertex has been the primary focus of this study, since it has 
the desirable properties of making the solutions approximately gauge-invariant and also 
multiplicatively renormalizable up to next-to-leading log order. Solutions have been obtained 
for comparison purposes, using the minimal Ball-Chiu vertex and using the bare vertex 
Ansdtze (with Zi = Z2). For the Ball-Chiu vertex, couplings in the range from = 0.1 
to 0.6 were used, while couplings up to oq = 1.75 were used with both the bare and the 
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Curtis-Pennington vertices. Various renormalization points and renormalized masses were 
studied. 

The subtractive renormalization procedure is straightforward to implement. The solu- 
tions are stable and the renormalized quantities become independent of regularization as 
the regularization is removed, which is as expected. For example, the mass function M(p^) 
and momentum-dependent renormalization Z{p'^) are unchanged to within the numerical 
accuracy of the computation as the integration cutoff is increased by many orders of magni- 
tude. For the range of couplings considered, the values of the renormalization constant 
Zi(/i,A) = Z2(/i,A) are never very far from 1 and vary relatively weakly with the choice 
of renormalized mass and cutoff as expected in Landau gauge. The mass renormalization 
constant A) converges with increasing A because the mass function M(p^) falls to 

zero sufficiently rapidly at large p^. This is a purely nonperturbative result and is in sharp 
contrast to the perturbative case where this constant diverges. For subcritical couplings, 
using the Curtis-Pennington vertex, we also find that the mass renormalization ^^(yU, A) 
scales approximately as A) oc (/i^/A^)*^2~'>cp(ao))^ where e.g. 7cp(0.5) = 0.358. 

Extensions of this work to include other gauges, other regularization schemes (e.g., soft 
cut-offs), and vertices of the Bashir- Pennington type are underway p4 . 



ACKNOWLEDGMENTS 



This work was partially supported by the Australian Research Council, by the U.S. 
Department of Energy through Contract No. DE-FG05-86ER40273, and by the Florida 
State University Supercomputer Computations Research Institute which is partially funded 
by the Department of Energy through Contract No. DE-FC05-85ER250000. This research 
was also partly supported by grants of supercomputer time from the U.S. National Energy 
Research Supercomputer Center and the Australian National University Supercomputer 
Facility. 



10 



REFERENCES 



[1] Johnson, N. Baker, and N. Willey, Phys. Rev. 136, Bllll (1964); ibid. 163, 1699 
(1967). 

[2] P. I. Fomin, V. P. Gusynin, V. A. Miransky and Yu. A. Sitenko, Riv. Nuovo Cim. 6, 1 

(1983). 

[3] J. E. Mandula, Phys. Lett. 67B, 175 (1977); A. Hey, D. Horn and J. E. Mandula, Phys. 
Lett. SOB, 90 (1978); J. Finger, D. Horn and J. E. Mandula, Phys. Rev. D 20, 3253 
(1979); J. Finger and J. E. Mandula, Nucl. Phys. B199, 168 (1982). 

[4] P. I. Fomin and V. A. Miransky, Phys. Lett. 64B, 166 (1976); P. I. Fomin, V. P. Gusynin 
and V. A. Miransky, Phys. Lett. 78B, 136 (1978); zbid., 91B, 421 (1980); 

[5] V. A. Miransky, Sov. Phys. JETP 61 (5), 905 (1985). 

[6] V. A. Miransky, Phys. Lett. 165B, 401 (1985); V. A. Miransky, Nuovo Cimento 90A, 
149 (1985). 

[7] A. Cohen and H. Georgi, Nucl. Phys. B314, 7 (9189); W. Bardeen, G. N. Leung, and 

S. T. Love, ibid. 323, 493 (1989). 
[8] U. Mahanta, Phys. Lett. 225B, 181 (1989). 

[9] H. Pagels, Phys. Rev. Lett. 28, 1482 (1972); T. Maskawa and H. Nakajima, 
Prog. Theor. Phys. 52, 1326 (1974); R. Fukuda and T. Kugo, Nucl. Phys. B117, 250 

(1976). 

[10] D. Atkinson and P. W. Johnson, Phys. Rev. D35, 1943 (1987); ibid. D37, 2290. 

[11] G. D. Roberts and A. G. Wilhams, Dyson-Schwinger Equations and their Application 

to Hadronic Physics, in Progress in Particle and Nuclear Physics, Vol. 33 (Pergamon 

Press, Oxford, 1994), p. 477. 
[12] P. E. L. Rakow, Nucl. Phys. B356, 27 (1991). 

[13] J. Oliensis and P. W. Johnson, Phys. Rev. D42, 656 (1990); K.-I. Kondo and 
H. Nakatani, Nucl. Phys. B351, 236 (1991); K.-I. Kondo, Nucl. Phys. B351, 259 (1991). 
[14] G-G. Wick, Phys. Rev. 96, 1124 (1954). 

[15] R. Haag and Th. A. Maris, Phys. Rev. 132, 2325 (1963); Th. A. J. Maris, V. E. Her- 
scovitz and D. Jacob, Phys. Rev. Lett. 12, 313 (1964); T. Nonoyama and M. Tanabashi, 
Prog. Theor. Phys. 81, 209 (1989). 

[16] J. G. Ward, Phys. Rev. 78, 124 (1950); Y. Takahashi, Nuovo Gimento 6, 370 (1957). 

[17] K.-I. Kondo, Y. Kikukawa and H. Mino, Phys. Lett. 220B, 270 (1989). 

[18] D. Atkinson, P. W. Johnson and K. Stam, Phys. Lett. 201B, 105 (1988). 

[19] J. S. Ball and T. W. Ghiu, Phys. Rev. D 22, 2542 (1980); ibid., 2550 (1980). 

[20] J. E. King, Phys. Rev. D 27, 1821 (1983). 

[21] P. Rembiesa, Phys. Rev. D41, 2009 (1990). 

[22] B. Haeri, Phys. Rev. D 43, 2701 (1991). 

[23] D. G. Gurtis and M. R. Pennington, Phys. Rev. D42, 4165 (1990). 

[24] D. G. Gurtis and M. R. Pennington, Phys. Rev. D44, 536 (1991). 

[25] D. C. Curtis and M. R. Pennington, Phys. Rev. D46, 2663 (1992). 

[26] D. C. Curtis and M. R. Pennington, Phys. Rev. D48, 4933 (1993). 

[27] Z. Dong, H. Munczek, and C. D. Roberts, Phys. Lett. 333B, 536 (1994). 

[28] A. Bashir and M. R. Pennington, Gauge Independent Chiral Symmetry Breaking in 

Quenched QED, University of Durham preprint DTP-94/48. 
[29] J. D. Bjorken and S. Drell, Relativistic Quantum Fields (McGraw-Hill, New York, 1965). 



11 



[30] F. Mandl and G. Shaw, Quantum Field Theory, (Wiley-Interscience, New York, 1984). 
[31] C. D. Roberts, Schwinger Dyson Equations: Dynamical chiral symmetry breaking and 

Confinement, in QCD Vacuum Structure, edited by H. M. Fried and B. Miiller (World 

Scientific, Singapore, 1993). 
[32] D. Atkinson, J. C. R. Bloch, V. R Gusynin, M. R. Pennington, and M. Reenders, Phys. 

Lett. 329B, 117 (1994) 
[33] C. Itzykson and J. B. Zuber, Quantum Field Theory, (McGraw-Hill, New York, 1980). 
[34] F.T. Hawes and A.G. Williams, in preparation. 



12 



FIGURES 



FIG. 1. Family of DSE solutions with the Curtis-Pennington vertex and subcritical couplings 
ao = 0.1 to 0.9. The renormalization point is /i^ = 100, and renormalized mass is m(/i) = 10 for 
all cases shown. Ultraviolet cutoffs are 1 X 10^2 foj, all 

cases, (a) the finite renormalizations A{p^); 

(b) the mass functions M{p'^). 



FIG. 2. Comparison of the solutions to the DSE with the three different vertices. All solutions 
have the subcritical coupling = 0.5, renormalization point /U^ = 100, and renormalized mass 
m(/u) = 10. Ultraviolet cutoffs are 1 x 10"*^^ for all cases, (a) the finite renormalizations A{p'^); (b) 
the scalar self-energy B{p'^). 



FIG. 3. Asymptotic scaling of the mass function for the three vertex Ansdtze at the subcritical 
coupling cto = 0.5; the functions plotted are M{p'^) x {p'^)^~'^ , where 7' = .3614329 is the correct 
anomalous dimension for the bare-vertex solution at this coupling. The anomalous dimensions 
extracted for the Ball-Chiu and Curtis-Pennington vertices respectively are 0.35836 and 0.35843. 



FIG. 4. The mass renormalization Zm as a function of /x^/A^, for the subcritical coupling 
ctQ = 0.5. Bare masses for all solutions were chosen to give mass m{p'^) = 10 at p"^ = 100. Boxes, 
□ , connected by solid lines, are results with cutoff A^ = 1 x 10^^; pluses, +, connected by dashed 
lines, are results with A^ = 1 x 10^^. Data are given in Table 0. 



FIG. 5. The scalar self-energy Bijp'), for supercritical coupling ao = 1-00 for the Cur- 
tis-Pennington vertex. Renormalization point is /i^ = 100, with renormalized mass m(/i) = 10. 
The cutoff is stepped by factors of 100 over eight orders of magnitude. Note that the mass curves all 
coincide. Above p^ = 2x 10^^, the mass curves are negative, but the maximum negative excursion 
is -2.4 X 10-"^. 



FIG. 6. Mass functions for the supercritical coupling oq = 1.15 with both the bare and Cur- 
tis-Pennington vertices. The renormalization point is /i'^ = 1 x 10*^, and m(^) = 400. 



FIG. 7. Zero-crossings in the self-energy BijP'). The case shown has qq = 1.75, /i^ = 1 x 10^, 
and m(^) = 500. Note that there are two zero-crossings, and that the function approaches zero 
from above as p^ 00. Note also that the negative peak is typically small compared with Bifi), 
m(/u), and p^. 
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TABLES 



TABLE I. Renormalization constants Z2(j[x, A), bare masses mo (A), and mass renormalization 



Zjn(iJ,,A) for various subcritical couplings for the Curtis- Pennington vertex, all with fj? = 100, 
m{n) = 10. Cutoffs A2 are 1 x lO^^ a,nd 1 x lO^^. 





A2 




mnfA) 


Z'm ill. A) 


0.1 


1 X 10^2 


0.99822798 


5.74672 


0.574672 




1 X 10^^ 


0.99982280 


4.0972 


0.40972 


0.2 


1 X 10^2 


0.99931590 


3.20125 


0.320125 




1 X 10^^ 


0.99931590 


1.59519 


0.159519 


0.3 


1 X 10^^ 


0.99851308 


1.7152 


0.17152 




1 X 10^^ 


0.99851308 


0.58245 


5.8245 X 10-2 


0.4 


1 X 10^2 


0.99744429 


0.874486 


8.74486 X 10-2 




1 X 10^^ 


0.99744428 


0.195803 


1.95803 X 10-2 


u.o 


± X lU 


n QQRl ^fiO^ 
u.yyoioozo 


n 41 7RQ9 


A 1 7fiQ9 V 1 n~2 




1 X lO^S 


0.99613623 


0.0590025 


5.90025 X 10-3 


0.6 


1 X 10^2 


0.99461286 


0.18240 


1.8240 X 10-2 




1 X 10^^ 


0.99461286 


0.015271 


1.5271 X 10-3 


0.7 


1 X 10^2 


0.99289572 


0.069783 


6.9783 X 10-3 




1 X 10^^ 


0.99289572 


3.14581 X 10-3 


3.14581 X 10-^ 


0.8 


1 X 10^2 


0.99100430 


0.0214029 


2.14029 X 10-3 




1 X 10^^ 


0.99100429 


4.37117 X 10"^ 


4.37117 X 10-^ 


0.9 


1 X 10^2 


0.98895625 


4.01753 X 10-3 


4.01753 X 10-3 




1 X 10^^ 


0.98895626 


2.31378 X 10-^ 


2.31378 X 10-^ 



TABLE II. Renormalization constants Zi and Zj^ for several different ratios j all solutions 
are for the subcritical coupling ckq = 0-5 for the Curtis-Pennington vertex, with renormalized masses 



chosen as m(/Lt = 


10) = 10. 










A2 = 1 X 10^2 




A2 = 1 X 10^^ 






Z2 




Z2 




1 X 10-1° 


0.99613623 


0.04177 


0.999999999439219 


0.032364 


1 X 10-9 


0.99914330 


0.04836 


0.999999999970789 


0.044846 


1 X 10-^ 


0.99993219 


0.06319 


0.999999999998478 


0.06214 


1 X 10-^ 


0.99999612 


0.08640 


0.999999999999921 


0.08611 


1 X 10-^ 


0.99999979 


0.11940 


0.999999999999995892 


0.11932 


1 X 10-5 


0.99999998 


0.16536 


0.999999999999999778 


0.16534 


1 X 10-^ 


0.9999999994 


0.22912 


1.000000000000000 


0.22911 


1 X 10-3 


0.99999999997 


0.31749 
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TABLE III. Renormalization constants Z2(/i, A), bare masses mo(/U,A), and mass renormal- 
izations Zm{^,K) for the supercritical coupling uq = 1.15 for the Curtis-Pennington vertex. All 
solutions are renormalized with fi^ = 1 x 10^, m{ii) = 400; the cutoff is stepped from 1 x 10^ to 

1 X 10l2_ 



A^ Z2 mo(A) Zrn 

1 X 10^ 1 0.999913461 230.581 0.57645 

1 X 10^ .1 0.999848292 53.5732 0.13393 

1 X 10^° .01 0.999846784 4.44299 1.1107 x 10"^ 

1 X 10^^ .001 0.999846902 -3.93251 -9.8313 x lO'^ 

1 X 10^2 _0001 0.999846922 -2.84691 -7.1173 x 10-^ 



15 



This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9410286vl 



This figure "fig2-l.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9410286vl 



o 

N 
• 1—1 

13 



o 



1.005 



1.000 



0.995 



0.990 




10^ 



10 



10 




p (spacelike) 



This figure "figl-2.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9410286vl 



This figure "fig2-2.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9410286vl 



< 1.005 

o 

^ 1.000 

• 1—1 

B 

g 0.995 

0.990 a) 
10 



\ 



10 



10^ 



10 



10 



o 

00 



10 
8 
6 
4 
2 




10" 









\^ bare 7^ 




\ minimal Ball-Chiu 




\ Curtis-Pennington 


renorm. point 




\i = 100 




m(|Li) = 10 


1 < 1 < 1 < 1 < 1 < 1 < 1 ,1 < 1 



10 



10^ 



10 



10 



p (spacelike) 



This figure "figl-3.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9410286vl 



Oh 



I 



Oh 



.54 



.52 - 



.50 - 



.48 - 



.46 - 



55 .44 - 



.42 - 



.40 



o 

00 



T ' — r 



T « — r 



T « — r 



T • — I « — r 




bare 7^ 

minimal Ball-Chiu 
Curtis-Pennington 



J I I I I I 



J I L 



J I L 



10' 



10 



10 



10 



12 



p (spacelike) 



This figure "figl-4.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9410286vl 



This figure "figl-5.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9410286vl 




10"^ 10° 10^ 10 



p (spacelike) 



1400 



1200 



1000 - 



800 - 



600 - 



400 - 



200 - 



-200 



-i-| 1-| i-i 1-| i-i i-i 1-| i-i 1-| 1-| i-i 1-| 1-| i-i 1-| rx 



bare 7^ 

Curtis-Pennington 



renormalization point 
[I = 10\ m(|Ll) = 400 



i-- 




10 



10 







10" 



10 



10 



p (spacelike) 




p (spacelike) 



